Non-commutative geometry as a realization of varying speed of light cosmology 
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We examine the cosmological implications of space-time non-commutativity, discovering yet another 
realization of the varying speed of light model. Our starting point is the well-known fact that non- 
commutativity leads to deformed dispersion relations, relating energy and momentum, implying 
a frequency dependent speed of light. A Hot Big Bang Universe therefore experiences a higher 
speed of light as it gets hotter. We study the statistical physics of this "deformed radiation", 
recovering standard results at low temperatures, but a number of novelties at high temperatures: 
a deformed Planck's spectrum, a temperature dependent equation of state w = p/ p (ranging from 
1/3 to infinity), a new Stephan-Boltzmann law, and a new entropy relation. These new photon 
properties closely mimic those of phonons in crystals, hardly a surprising analogy. They combine 
to solve the horizon and flatness problems, explaining also the large entropy of the Universe. We 
also show how one would find a direct imprint of non-commutativity in the spectrum of a cosmic 
graviton background, should it ever be detected. 



I. STATEMENT OF PURPOSE 

Inflation and string theory are widely perceived as the 
leading schools of thought in cosmology and fundamen- 
tal physics, respectively. And yet finding common ground 
between them, say by deriving inflation from a concrete 
string theory model, has so far remained elusive. Such an 
enterprise is far from frivolous, since cosmology is proba- 
bly the most realistic laboratory for testing string theory. 
For this reason it is disappointing that after all these 
years, and following a number of radical revolutions in 
both fields, the two theories have still failed to condense 
into a single construction. 

Could it be that somehow we got things wrong on one 
side of the story? The inspirational origins of infiation 
are founded on a number of problems of Big Bang cos- 
mology, namely the horizon, fiatness, homogeneity, and 
entropy problems [1-4]. A recently proposed alternative 
explanation to these problems is the varying speed of 
light (VSL) theory [5-7], the idea that fight travelled 
faster in the early Universe. Currently VSL is far less 
developed than inflation, with full model building work 
hardly started. But could VSL be more amenable to a 
direct connection with fundamental physics [8,9], more 
specifically with string theory? 

In this paper we establish a connection between non- 
commutative geometry and VSL. Over the last few years 
non-commutativity of space-time coordinates has become 
part and parcel of any attempt to quantize gravity. Simi- 
lar behaviour was already spelled out in Appendix I of [6] 
(see also [7]), where it was shown that under VSL the par- 
tial derivatives dt and dx do not commute. Furthermore 
non-commutative geometries arc known to lead to de- 
formed dispersion relations [20] , which imply a frequency- 
dependent speed of light. 

A deep connection is up for grabs, and such is the 
purpose of this paper. 



II. COSMOLOGY AND NON-COMMUTATIVITY 

One of the pillars of the standard Big Bang (SBB) 
is the theory of General Relativity. The problems of 
the SBB model suggest that classical General Relativ- 
ity is inconsistent in the very early universe where quan- 
tum corrections to gravity become significant. Interest- 
ingly, most attempts toward a theory of quantum gravity 
such as string/M-theory realize that space-time itself is 
non-commutative. In particular, non-commutativity is 
concretely realized in a non-perturbative formulation of 
M-theory, (M)atrix theory, where the D-0 brane collec- 
tive coordinates are in general non-commuting; this is 
interpreted as a non-commuting space-time. In the infi- 
nite momentum frame (IMF), non-commuting space-time 
emerges from (M)atrix theory [11]. Also, in the open 
string sector of string theory, non-commutative geome- 
try arises when the NS-B field is turned on; the open 
string end points become "polarized" and hence non- 
commuting [13,12] . In the low energy limit of open string 
theory in the presence of a B field, these dipoles are re- 
alized as non-commutative solitons in Non-commutative 
Yang-Mills (NCYM) theory where they exhibit a VSL 
behavior [14]. 

The coordinate operators of a non-commuting space 
obey the following commutation relation. 



(1) 



where O^'' is in general constant, antisymmetric, and 
with dimension [length]"^. There exists a general pre- 
scription for studying quantum field theories in non- 
commutative spaces with flat backgrounds. By simply 
taking a classical field theory defined on a commuting 
space and replacing the product operation by a * prod- 
uct, one obtains a non-commutative field theory. The * 
product is defined from 



(2) 
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= + '-e^^^d^^d,^ + o{Q^) (3) 

Notice that the * product is non-commutative. Also non- 
commutative field theories are necessarily non-local since 
the * product yields an infinite number of derivatives in 
the Lagrangian. This non-locality in space-time, as we 
will see in the sections that follow, gives rise to a VSL 
due to a modification in the dispersion relation of the 
quantum fields in a non-commutative space-time. 

Although non-commutativity is concretely realized in 
string theory, space-space non-commutativity is better 
understood than its space-time counterpart. In both 
cases, however, non-commutativity is best understood 
for constant B field (theta parameter); while more com- 
plicated versions of non-commutativity in string the- 
ory generalized to non-constant B-field and also its 
gravitational analogue is currently under investigation 
[16,15,12]. For VSL we arc concerned about space-time 
non-commutativity with non-constant B field (Q). 

In light of the string-theoretic realization of non- 
commutativity, we shall implement a more general alge- 
bra, where the parameter is not constant, but a linear 
function of x, a quantum deformed Poincare' group. This 
extension is important because it precisely yields a VSL 
in the photons by deforming their dispersion relations. In 
equation (1) relativistic invariance remains classical and 
the translations 

^ ^f_i ~ ~t~ (4) 

preserving the algebra of (1) are commutative. However, 
if depends on x^jin particular if we assume the linear 
x^ dependence, i.e. 

= ic^j^XA- (5) 

then, the translations 

x'i^^x^+ (6) 

are not commutative. For particular choices of the 
structure constants c^^ this type of non-commutative 
structure falls under the category of the quantum- 
deformed Poincare' group, with non-commutative group 
parameters^. In this work we consider a particular 
type of quantum Poincare' group; the k deformation of 
space-time symmetries with non-commutativity of space- 
time coordinates described by the following solvable Lie- 
algebra relations: 

1 

[xo,Xi] = -Xi, [xi,Xj] = 0. (7) 



^For a more descriptive classification of the quantum 
Poincare' groups see [10] 



Thus spatial translations still commute; 
non-commutativity affects only operations involving the 
time coordinate. 

This type of non-commutativity has yet to be realized 
in string/M theory, however we expect that it will be 
better understood in the future, especially in the con- 
text of the non-commutative version of AdS/CFT corre- 
spondence and in formulations where the B field (0) is 
non-constant [14,15]. 

In the case of the deformed Poincare algebra, the spa- 
tial momenta remain commutative 

[P.,P,]=0 (8) 

and the rotation part of the Lorentz sector is left un- 
changed. For this reason we assume there is a scale Ajvc 
in the Early universe such that 



where k is 0(1); ie. non-commutativity is strong. As a 
result it is safe to conjecture that since 

(1) The K deformation does not act on the rotations 
and the canonical momenta. 

(2) The deformation is strong, k ~ 0(1) , 

thus the non-commutativity can be described as a gas 
of massless non-commutative radiation propagating in a 
conformally flat Friedmann-Lemaitre-Robertson- Walker 
space-time. 

The radiation is in thermal equilibrium and is maxi- 
mally correlated. This is a statement that the space-time 
is a Non-commutative FRW (NCFRW) with a gas of pho- 
tons. The modified dispersion relation of the photons 
captures the non-commutativity of the space-time. The 
FRW space-time is endowed with fiat spatial sections, 
and is given by the line element 

ds2 = a^(j^){-dff + dx^ + dy^ + dz^). (10) 

where r] is the conformal time and {x,y,z) are spatial 
comoving coordinates. This metric is conformally flat 
and is locally equivalent to the Minkowski metric r^^^. 

Since our space-time is non-commutative, the notion of 
a continuous space-time manifold is lost. In fact the sep- 
aration between energy-momentum and curvature, which 
is manifest in the Einstein field equations is not as trivial. 
However, we can transcend this difficulty by deforming 
the dispersion relations of the photons and study their 
dynamics in a commuting space-time; this yields the na- 
ture of how matter behaves in a quantized space-time. 
So long as the photons are in equilibrium, our approxi- 
mation is valid. 

Therefore, our starting point for cosmology is a de- 
formed Poincare group whose action does not affect the 
conformal factor of the metric. A similar approach was 
taken in the context of non-commutativity and inflation 
[17]. This is analogous to the standard hot big bang 
scenario where our light degrees of freedom consist of a 
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gas of non-commutative radiation in thermal equilibrium. 
This radiation will exhibit a modified dispersion relation 
which alters its equation of state. 



III. THE THERMODYNAMICS OF 
NON-COMMUTATIVE RADIATION 

In the rest of the paper we will study the cosmologi- 
cal consequences of modified dispersion relations for pho- 
tons. However, it is instructive to provide a discussion of 
how modified dispersion relations arise in the free field 
regime of the deformed n Poincaire group. In what fol- 
lows, we will provide one of a few examples of a modified 
frequency dependent dispersion relation. 

In [19] the authors were able to define at the perturba- 
tive level a free field theory for massless bosons by iden- 
tifying the propagator and constructing the intertwiners 
between the representations of the Poincare group act- 
ing on fields and states. These intertwiners are simply 
the wavefunctions of particles of definite spin. One can 
define these intertwiners through the appropriate wave 
equations describing free fields. 

Therefore, the dispersion relations can be obtained 
from the invariant wave operator on the K-deformed 
Minkowski space -^r- , which is expressed in momen- 
tum space as 



CI 



bcp 



4k2 



where 



^bcp ^ (p2)g-Po/« _ (2Ksinh(J^^ 



(11) 



(12) 



is the first Casimir of the K-deformed Poincaire algebra. 

It follows that the spectrum of the modified massless 
wave operator contains the deformed massless mode 



Cl^PcP = 

which leads to the following dispersion relation 



1 K 



\2Ksinh\—^ = 



(13) 



(14) 



The wave equation leading to the above dispersion rela- 
tions is non-local in time but can be reexpressed so that 
the dispersion relations corresponds to an operator of 
second order in time derivatives and non-local in space. 



a;2 = 



k 

Klog{l H — ) 



(15) 



Here f{E) gives rise to a frequency dependent speed of 
light, whereas c is a possibly space-time dependent, but 
frequency iiiciepenclcnt, speed. We thus factorize two dif- 
ferent VSL effects previously considered in the literature - 
a frequency dependent c (as studied in [20]), and a space- 
time dependent c (as examined in [5-7]). Most likely, the 
final theory will contain both. Naturally the most gen- 
eral deformation need not factorize these two effects, or 
even produce a quadratic invariant as in Eqn. (16), but 
we shall retain this assumption. 

We shall consider two proposals for /, previously con- 
sidered in the literature. Amelino-Camelia and collabo- 
rators (eg. [21]) have considered 



f = l + \E 



(17) 



associated with the K-Poincare group, whereas following 
the quantum deformation of the Poincare group as stud- 
ied by Majid one finds [20]^: 



/ = 



2XE 



-2XE 



(18) 



Here A is the deformation parameter and we have set 
h = 1. The last deformation has also been proposed by 
Kowalsky-Glikman. We will also consider a generaliza- 
tion of (17), of the form 



/ = (1 + A^) ' 



(19) 



in particular for 2/3 < 7 < 1. The case 7 > 1 will be 
studied in detail elsewhere [32]: it leads to a realization 
of inflation, but not VSL. We shall call these models 1, 
2, and 3, respectively. 

At once we find two possible definitions for the speed 
of light. One may use the definition proposed in [20]: 



cf 



dE 



dp 1 - ^ 



(20) 



(where ' denotes a derivative with respect to E) for which 
model 1 gives: 



c = c{l + XEf 
whereas model 2 gives 



■ ce 



2XE 



Alternatively we may define the speed of light as 

c = — =cf . 

P 



(21) 



(22) 



(23) 



Both definitions, c and c, play a role in what follows. 



Let us now consider the following generalization of 15 
for massless particles, for which: 



E' 



2 2 j-2 







(16) 



■^Wc have applied A — 2A to the convention used in [20] so 
that the 2 dispersion relations to be considered agree to first 
order 
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A. The partition function and phase space densities 

We now proceed to examine the statistical physics of 
this "deformed" radiation. We first note that the form 
of the partition function 



(24) 



does not depend on the dispersion relations; indeed it 
amounts to a definition of temperature based upon the 
fact that the thermal reservoir has a much larger number 
of states than the system under examination. Hence, for 
a boson gas, the average number density of states with a 
given momentum still satisfies the Bose-Einstein result: 



n(p) 



1 



g/3S(p) _ I 



(25) 



since this formula depends only on the partition func- 
tion and the rules of counting associated with bosons. A 

similar undeformcd expression applies to fcrmions. 

The density of momentum states for point particles is 
again unchanged, as it simply reflects the use of periodic 
boundary conditions (p = (2Trh/L)n where n is a triplet 
of numbers, and L is the side of a given cubic volume). 
All that changes is the relation between E and p and 
therefore the density of states Q{E) with a given energy 
E per unit of volume. We find: 



n{E) = 



E^ 



1-^1. (26) 



simply from Q{E)dE = Q{p)dp. 



B. The deformed Planck spectrum and the graviton 
background 

From p{E) = n{p{E))fl{E) we thus obtain the de- 
formed Planck's spectrum: 



p{E) = 



E^ 



1- 



f'E 



(27) 



One can easily show that only deformations of the form: 

f = [1 + {XE)-^]-^/^ (28) 

leave the Planck spectrum unchanged. These are not 
very physical in a cooling Universe, where the speed of 
light would initially be a constant, then start to drop, 
with c oc l/o, until nowadays. 

All other dispersion relations affect Planck's law. This 
does not conflict with experiment if / ~ 1 for T <C 1/ A. 
However even then there may be an observational im- 
print. Gravitons decouple around Planck time, and 
should constitute a thermal background similar to the 
cosmic microwave background (with a temperature of the 



same order). However their spectrum will be a deformed 
Planck spectrum, since it mimics the spectrum they had 
at the time they decoupled. Hence if and when a gravi- 
ton background is discovered its spectrum will supply a 
direct measurement of the dispersion relations. 

More specifically model 1 leads to the deformed spec- 
trum: 



p{E) = 



1 



1 



TT-^K^C^ 6/3^ -1(1 + XEY 



(29) 



For T <C 1/A distortions to this spectrum are negligible. 
For high temperatures, though, the peak in p{E) is not 
at T, but at = 1/A for all temperatures. In other 
words the color temperature saturates at T^ax ~ 1/A. 
The general form of the spectrum in this regime is 



p{E) = 



E^T 



-K^^{l+\EY 



Model two leads to: 
p{E) = ^ 



E/X^ 



-i\E 



sinh'^iXE) 



which for T > 1/A leads to 

T 



p{E) 



XH^ 



sinh^{XE) 



(30) 



(31) 



(32) 



In both cases the spectrum is super-black (as opposed to 
gray); that is, the overall amplitude of the spectrum with 
respect to black at T « 1/A is enhanced by a factor of 
TX. In Fig. 1 we plot spectra for model 1 and 2, as well 
as the undeformcd Planck spectrum. 

Model 3 is unique in that it never saturates the color 
temperature, i.e. for all temperatures hotter radiation 
means a peak in the photons distribution at higher ener- 
gies. It's thermal spectrum is given by 



p{E) 



1 E"" (1 + (1-7)(A£;)T') 
7r2n^c3 - 1 (1 + (A^)7)4 



(33) 



For 7 < 2/3 the peak of p{E) scales like T. For 2/3 < 
7 < 1 the peak becomes very wide and covers energies 
from E K. X~^ to E p^i T. We illustrate this feature in 
Fig. 2. 



C. The energy and entropy densities 

Integrating p{E) leads to a modified Stephan- 
Boltzmann law. For low temperatures p oc T"* as usual, 
but for AT 1 we find that p oc T. This dependence 
can be expressed in terms of the function 



n{E) 



dlogp 



(34) 



rflogT 

which we plot in Fig. 3. At low temperatures the pro- 
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1 Planck aL '1'= 1 /A 




FIG. 1. Deformed thermal spectra at high temperature 
(TA 3> 1), divided by (we have assumed A = 1). For refer- 
ence we plotted a Planck spectrum at temperature T = 1/A. 
For model 1 and 2 the color temperature saturates, so that 
above T = 1/A the peak in p{E) does not shift to higher 
energies. 




FIG. 2. Deformed thermal spectra for model 3, with 
7 = 1/2, as the temperature is increased. We see that for 
this model at all temperatures, hotter radiation implies higher 
mean energies for photons, and a shift to the higher energies 
of p{E). 
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FIG. 3. The transition from p oc T* to p oc T at very high 
temperatures for non-commutative radiation. 

portionality constant is the usual Stephan constant. At 
high temperatures the proportionality constant is model 
dependent. In general, at high temperature: 

, T 



A3 



with 



a+ = 



1 

3^ 
1 

247r2 



(35) 



(36) 



(37) 



for models 1 and 2, respectively. 

For model 3 the situation is a bit different. The 
high energy Stephan-Boltzmann equation takes the form 
p oc with a power ( in the range 1 < C < 4. An 
approximate high temperature formula for C may be ob- 
tained by noting that in Eqn. 33 the last factor is approx- 
imately 1/(AT)3'>' over the most relevant parts of the inte- 
grand. Hence we may expect at high energies p oc T^-^^ 
that is C « 4 — 37. We have verified numerically that 
this analytical approximation works extremely well. In 
Fig. 4 we describe the transition from the low to the high 
temperature behaviour for various values of 7. 

The entropy density then follows from 

dTdp_ 
T dT 



(38) 



For low temperatures we recover the usual expression 

4 



s = -aT-" 
3 



whereas for high temperatures we now have 
s w 0-+ log AT 



(39) 



(40) 
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-2-10123 

log AT 

FIG. 4. The transition from p oc to p oc at very high 
temperatures for model 3. Wo see that for 7 < 1, the high 
temperature relation has a power in the range 1 < C < 4. In 
all cases hotter radiation translates into denser radiation. 
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log (pA*) 

FIG. 5. The equation of state for deformed radiation for 

model 1 and 2 (solid lines). The dashed lines represent fits 
of the form w = i + log(l + Dlog(pX'^)) with D « 4 and 
D fs 250 respectively, and the dotted line shows the slope 1. 



D. The equation of state 

The pressure formula (and so the equation of state) is 
also modified. As before it is given by a sum over all 
states 



n. 



dV 



but since now we have for each state 



(41) 



(42) 



we have 



1 ^ 1 [ p{E)dE 1 

^ / / 

and so the equation of state of radiation is modified. 

To find the modified equations of state we therefore 
have to compute the integral: 



1 f p{E,T)dE 



(44) 



Combined with p{T) this leads to a modified equation 

of state p = w{p)p. It is not difficult to guess its form. 
At low temperatures w « 1/3. At high temperatures we 
have 



P^^{A + \og{TX)) 



(45) 
(46) 



where A is order 1. So, for TA » 1, we have 

w{p) + log(TA) kB + log(pA^) (47) 

that is w grows logarithmically with p (and therefore with 
T) . This approximate argument is confirmed by a numer- 
ical integration as shown in Figure 5. 

Not siirprisingly one may violate all of Hawking's en- 
ergy conditions apart from the weak one (the energy 
density is still always positive). Indeed a sign of non- 
commutativity appears to be that it generates equations 
of state with \w\ > 1. The fact that the usually sensible 
energy conditions are badly violated should not deter us. 

Model 3 leads to a simpler equation of state. At 
low energies w = 1/3; at high energies the equation of 
state becomes again a constant. Indeed for this model 
{f'E)/f = J at E A^^ so that the denominator of 
Eqn. 44 becomes approximately a constant at high tem- 
peratures, over the peak of p{E). Hence we may expect: 



Wo 



= w{p 00) 



3(1-7) 



(48) 



A numerical integration reveals that this is only an ap- 
proximate formula, albeit with the correct qualitative be- 
haviour: in Fig. 6 we show the asymptotic equation of 
state as predicted by our analytical fit and by a numeri- 
cal integration. In [32] we shall explore the dramatically 
different range 7 > 1 , to find that it is possible to gener- 
ate an inflationary equation of state. 

Notice that the results in this section lead to a closed 
form expression for the entropy: 



P + P 



(49) 
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FIG. 6. The equation of state for T ^ for model 3 with 

7 < 1. The dashed hnc represents the analytical fit described 
in the text; the solid line the result of a numerical integration. 



with limiting cases (39) and (40). 



E. The phonon analogy 

We conclude this section by noting that to some ex- 
tent the thermodynamics of deformed radiation is very 
similar to that of phonons in crystals. Thermal phonons 
satisfy Bose-Einstein statistics; however, unlike photons, 
they are subjected to very complicated dispersion rela- 
tions. Hence the density of states for phonons Q.{E) (cf. 
Eqn. (26)) at high energies is very different from that of 
photons (see for instance Fig. 8-4, pp 217 of [22]), giving 
rise to a multitude of thermodynamical novelties simi- 
lar to the ones we have just derived for deformed radi- 
ation. At low energies the phonons cannot sec the dis- 
crete structure of the crystal and behave like ordinary 
photons. At high energies, on the contrary, they become 
highly sensitive to crystal properties and exhibit exotic 
behaviour. The borderline between these two regimes, for 
thermal phonons, is determined by the Debye tempera- 
ture, which is therefore a solid state physics counterpart 
to the Planck temperature. 

The analogy we have spelled out is far from surpris- 
ing. After all space-time non-commutativity is nothing 
but a method for quantising space-time. In some sense 
this means introducing a discrete space-time structure 
not dissimilar to that of a crystal (another interpretation 
is space-time uncertainty [29]). As photon frequencies 
get higher and higher (or their temperature approaches 
the Planck temperature) their dispersion relations, and 
all derived thermodynamical properties, start being sen- 
sitive to the discreteness of the space-time structure sup- 



porting them. 



IV. THE TRANS-PLANCKIAN COSMOLOGICAL 
EVOLUTION 

We now proceed to integrate the Friedmann equations 
with the w = w{p) peculiar to non-commutative radia- 
tion. We first consider the case of zero spatial curvature 
K = 0, leaving for Section VI the cases K = ±1. The 
equations may be written^: 



1 



- = --p{i + Mp)) 

a 6 



(50) 

(51) 



where a is the scale factor and dots represent derivatives 
with respect to proper time. Instead of the second equa- 
tion one may use the conservation equation 



p + 3-{l + w{p))p = 



(52) 



From the equations above it is also possible to find an 
equation for the speed of sound Cg, defined by: 



2 _dp _p 
Sp p 



(53) 



Computing w and using equation (52) leads to [24]: 



w 



-3-{l+w){ci-w) 



and so 



w + p 



dw 
dp 



(54) 



(55) 



an equation which shall be of great relevance in a future 
publication, in which we discuss density fluctuations in 
these scenarios [23]. 

For model 3 it is easy to obtain an analytical solution 
in the asymptotic trans-Planckian regime: 



(56) 



A numerical integration (with a w{p) also obtained nu- 
merically) is presented in Fig. 7. We see that after 
Planck's time the imivcrsc expands as a normal radiation 
dominated universe, but inside the Planck epoch it ex- 
pands slower, the more so the higher the value of 7. This 
is because the decceleration of the universe is higher: 



aa 



1 2 



2 1 



(57) 



*In what follows we shall set 8itG = 1 



7 




FIG. 7. The cosmic evolution (with t and p is Planck units) 
for Model 3 with various values of 7. Model 1 can be in- 
ferred from 7 — > 1. The expansion is always decelerated, even 
though the universe expands slower inside the Planck epoch. 



Nowhere, in the transition between one regime and the 
other, do we have accelerated expansion. This is to be 
contrasted with the case 7 > 1, explored elsewhere [32]. 

Models 1 and 2 asymptotic trans-Planckian behaviour 
may be inferred from (47) and (52). For pX* ^ 1: 



expj^ 



(58) 



where C is a constant. Integrating (50) then leads to 



Ei 



2^ 



oc t 



(59) 



The qualitative behaviour is similar to the one described. 
A numerical intergation leads to a result equivalent to 
model 3 with 7 — > 1. 



V. A SOLUTION TO THE HORIZON PROBLEM 

Naively it may seem self-evident that the horizon prob- 
lem has been solved in this model. As the Universe gets 
hotter, the photons' average frequency and energy get 
higher, just like in the standard Big Bang model. How- 
ever in our theory, unlike in the standard model, this 
results in a larger "ambient" speed of light. Hence our 
model realises the usual VSL solution to the horizon 
problem [6], not by a direct time-dependence of c, but 
indirectly, via the chain linking time, temperature, aver- 
age photon frequency, and c. 

However, closer scrutiny shows that it is not that sim- 
ple; in fact not all deformations / lead to a solution to the 



horizon problem. On the one hand our analysis confirms 
that as t decreases below Planck's time, the temperature 
keeps increasing. This is obvious for model 3 but not for 
model 1 or 2. However we see that a keeps decreasing 
(cf. Fig 7 and Eqn (59)), and from Eqn. (58) and (35): 



T = 



(60) 



implying that the temperature T does diverge in the early 
Universe for all models considered. 

However, as T grows above , we find that for model 
1 and 2 (but not for model 3, with 7 < 1) the color 
temperature T^. (defining the peak of p{E)) saturates at 
Tc = A~^. Hence a hotter plasma does not contain more 
energetic photons at the peak of the distribution; it only 
contains more photons (proportionally to T^) at a peak 
located at the same energy. Hence the most abundant 
photons in pre-Planck times will experience a roughly 
constant speed of light. Given that we do not have ac- 
celerated expansion, at first it looks as if these two VSL 
scenarios do not actually solve the horizon problem. 

Fortunately a further subtlety comes into play. Recall 
that for T ^ we have: 



p{E) = 



1 E^T 



1-^) (61) 



up to £■ Ri T. For E > T the distribution is then domi- 
nated by the exponential cut off e~^^: 



p{E) 



1 



fE_ 
f 



(62) 



For model 2, p{E) falls off exponentially away from the 
peak at Tc ~ A~^, but for model 1 the fall-off is merely 
power-law. More concretely for model 1 there is a non- 
negligible density of photons with E a T, the fastest 
photons in the Universe, given by 



p{T)T : 



1 + {XTY 



(63) 



The conclusion is that even though the color temperature 
saturates for T 3> A~^, there is still a constant density (of 
the order of the Planck density) of photons with energy 
of order E pa T. These are the fastest photons in the 
Universe, experiencing a speed: 



c = c{l + XT) = c{l 



AV 



(X exp 



(64) 



or a similar expression for c. Fast photons at the tail at 
E = T ensure causal contact at a given time. It is there- 
fore essential that interactions exist between photons at 
E = T a,nd at E = X-^. 

A similar calculation for model 2, on the other hand, 
leads to: 



p{T)T « Te 



-2AT 



(65) 
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FIG. 8. A log-log plot of the comoving horizon distance 
r/i = c/d as a function of the expansion factor a. At late 
times, as c stabilizes and a acquires the standard a oc t^^^ 
behaviour, we recover oc a. Before Planck's time, however, 
the decreasing c (or c) forces rn to decrease, thereby solving 
the horizon problem. 



implying that the relevant photons become exponentially 
suppressed. Model 3 on the other hand does not suffer 
from any of these problems. Nonetheless the conclusion 
remains that not all deformations / load to a sohition to 
the horizon problem in this VSL scenario; but some do. 
Let us concentrate on model 1 and 3. 

The causal range at time t in model 1 is defined by 
the distance travelled by photons with E during one 
expansion time (the period over which expansion can be 
neglected). The comoving causal range is therefore: 



cH- 



c 
a 



(66) 



At late times (when TX <C 1), c stabilizes, and a oc t^/'^, 
leading to the growth of the comoving horizon ru <x. a. 
Hence the comoving region we can now see was once split 
into many disconnected regions. For TA 1 we find 
instead: 



C 

- cx 



p exp ( 
cx ^ 



c 

2a^ 



(67) 



a decreasing function of a: hence providing a solution to 
the horizon problem. 

In Fig. 8 we plotted r/j as computed from a numerical 
integration for model 1. We have plotted the horizon as 
defined by c and by c. In both cases we see the distinc- 
tive late time growth Vh oc a, but the early exponential 
decrease in r/,. We have plotted how a given comoving 
scale is initially in causal contact, then falls out of causal 



contact as the speed of light slows down, to reenter the 
Hubble radius later on. 

For model 3 the solution to the horizon problem is 
more straightforward. Since the color temperature never 
saturates, as radiation gets hotter so does the ambient 
speed of light; we find that 



c oc c cx (AT)i 



(68) 



without any extra assumptions. Hence, using p cx 
and (56) we find 



37-2 

rh=T— 



(69) 



so that for 7 > 2/3 we solve the horizon problem. 



VI. A SOLUTION TO THE FLATNESS PROBLEM 

The discussion in Section II leading to the establish- 
ment of a non-commutative Friedman metric, and equa- 
tions, breaks down if wc introduce spatial curvature 
K = ±.1. In that case our space-time is conformally 
related not to non-commutative Minkowski space-time, 
but to a non-commutative (pseudo-)sphere, namely the 
fuzzy three-sphere S^. For a derivation of see [18]. 

We know that such a fuzzy sphere has the usual prop- 
erties of ordinary spheres as long as the curvature radius 
is much larger than the Planck length. However as one 
tries to curve them beyond Planck curvature, their cur- 
vature radius saturates. Hence we avoid a singularity. 
This picture, however, will not provide a solution to the 
flatness problem. Effectively wc may model this effect 
by replacing the curvature K/a? by Kg^/a?', where g{a) 
describes a deformation of the curvature of the sphere as 
it approaches Planck scale. A suitable function is: 



9{a) = 

and we may check that: 

Kg' 



1 



1 + 



1 

A2 



(70) 



(71) 



for a ^ A. It is clear that this model would convert the 
curvature term into a cosmological constant term, leading 
to a pre-Planck deSitter phase which then decays into a 
curvature dominated phase. 

We therefore adopt a different approach: we model the 
fuzzy sphere's unusual curvature by means of a direct 
coupling between K and the matter density p. The idea 
is that the curvature of spatial surfaces depends upon 
the energy scale of the matter probing it. This creates 
an intertwining between matter and geometry, which is 
to be expected above Planck scale. Indeed one does not 
expect the usual division between matter and curvature 
to survive a quantum gravity epoch (see discussion in 
the second to last paragraph of Section II). Hence we 
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replace the curvature term K/a^ by Kg"^ / , where now 
g = g(p) describes a deformation of the curvature of the 
sphere as the matter filling up the Universe approaches 
Planck scale. We choose: 



g{p) = {l + pX'y 



(72) 



As it turns out this approach provides a direct real- 
ization of the solution to the flatness problem described 
in [6] (just replace the c in Friedmann equations by 
Ck = cg{p))- The Friedmann equations are now: 



3^ 



Kg'6 



,2„2 



which lead to the integrability condition: 

^ g P 



p + Z-{l + w{p))p = 



(73) 
(74) 

(75) 



(where ' means a derivative with respect to p) or instead: 
QKc^gg'' 



1 



-3-1 
a 



■ w{p))p 



(76) 



We find therefore a coupling between the curvature K 
and the matter density. In order for flatness to be stable, 
g' should be positive: then supercritical models {K = 
1) have matter removed from them, sub-critical models 
{K = —1) have matter put into them, thereby creating a 
flatness attractor. 

Let pc be the critical density of the Universe: 



Pc = 



SttCVo 



(77) 



that is, the mass density corresponding to = for 
a given value of a /a. Let us define e = — 1 with 
n = p/pc. A numerical integration of equations (73) 
and (75), expressing the result in terms of e(a), is plot- 
ted in Fig. 9. We see that at late times, when the effects 



of non-commutativity have switched off, 



(the flat- 



ness problem), but early on e decays very fast, producing 
a very flat Universe at the end of the Planck epoch. 

This result can be analytically understood by combin- 
ing Eqns. (73) and (75) as in [6] into a single equation 
for e: 



e = {l + e)e-{l + 3w) + 2-e 
a g 



Hence, for e <C 1, and for pA^ 1, we have 



(78) 



(79) 



showing that any exponent a > 1/2 leads to a solution 
to the flatness problem. 



a = 2, X=l 



o 



log a 

FIG. 9. A log- log plot of |e| versus a. We see that at late 
times, when the effects of non-commutativity have switched 
off e oc (the flatness problem), but early on e decays like a 
power of p. 



For illustration purposes in Fig. 9 we started the inte- 
gration at T « lO^Tp, and with e ~ 1. In our scenario, 
however, there is no starting time, subjected to "natural" 
initial conditions (we recall that in the Big Bang model 
one likes to impose e ~ 1 a.t t — tp). Our model plunges 
as deep as we like inside the Planck epoch, and as we 
have shown, we found that during this period flatness 
then becomes an attractor. Should this regime be valid 
all the way up to T — oo we may conclude that when 
Planck time is finally reached e equals precisely zero. 



VII. THE ENTROPY OF THE UNIVERSE 

The flatness problem is often rephrased as an entropy 
problem, both in its infiationary and VSL formulations. 
As explained in [1,27,6] adiabatic solutions to the flat- 
ness problem lead to a seemingly very flat Universe, but 
containing only one particle within each curvature ra- 
dius' volume. This problem afflicts Brans-Dicke based 
solutions to the flatness problem [26,27] 

All our work assumes thermal equilibrium, so it would 
seem that we have an adiabatic scenario - that is, no 
entropy is produced. However this is not the case: we 
have introduced source terms in the energy conservation 
equation (cf. Eqn. (76)), and this necessarily modifies 
the first law of thermodynamics, a fact which, as we now 
show, leads to entropy generation even without leaving 
thermal equilibrium. 

Taking the equilibrium expression 

Td{sa^) = d{pa^) + pda^ = a^dp +{p + p)da^ (80) 
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one may read off from (75) that 



that is: 



S + 3-S = 

a 



GKc'^g . 



(81) 



(82) 



We see that ii K = entropy is conserved, a fact which 
may be checked directly from Eqns. (40) and (60). How- 
ever, if K = —1 and a > 0, then entropy is produced 
even though we have only used equilibrium thermody- 
namics. The case K = 1 and a > seems to lead to 
entropy reduction, and thus it may be argued that it is 
inconsistent with the second law of thermodynamics (see 
[28] for such an argument in a difFcrcnt context). 

We may now define the entropy (or number of parti- 
cles) inside a volume with the curvature radius as 



CTK = s 



and from (82) we have: 



6a 



(83) 



(84) 



As long as K = —1 and a > we see that ax increases 
unboundedly, under general conditions, before Planck 
time. After Planck time it stabilises to a constant. As 
in the discussion at the end of Section VI we therefore 
conclude that the natural state at the end of the Planck 
epoch is (Jk = oo, explaining the current bound that gk 
must be bigger that 10^*. 

The entropy problem is sometimes referred to as the 
need to explain why the Universe is so big, or why it 
contains so much entropy. Usually one needs to invoke 
an entropy production mechanism, such as reheating at 
the end of inflation, in order to solve this problem. In 
our scenario the explanation is related to the fact that 
in the non-commutative phase of the Universe we break 
Poincarc invariance. We know that energy conserva- 
tion follows from time-translation invariance - a property 
which is now deformed. Hence we expect violations of en- 
ergy conservation, and these lead to entropy production 
even if we remain in thermal equilibrium. 



VIII. CONCLUSIONS 

There are many avenues which lead to the possibility 
that space-time at the Planck scale is non-commutative. 
In the context of the Early Universe, VSL is able to 
solve the outstanding problems of SBB. In this paper 
we demonstrate that non-commutativity, in the context 
of the early universe, is equivalent to VSL. We provide a 
concrete cosmological model of a non-commutative early 



universe scenario which exhibits VSL phenomena. In this 
model the speed of light does not depend explicitly on 
time; rather it depends on the photon's frequency, and 
so as the Universe gets hotter the "ambient" speed of 
light increases. We showed that in this model the hori- 
zon and flatness problems are resolved. 

Perhaps the most important result of all, however, is 
the discovery that our solution to the horizon and flatness 
problems is related to the emergence of Lorcntz invari- 
ance in the Universe. We found that photons in non- 
commutative space-times are very similar to phonons in 
crystals, with Planck's temperature playing the role of 
the Debye temperature (see Section HIE). At low en- 
ergies the photons perceive a Lorentz invariant contin- 
uum, and have a constant speed. This corresponds to 
the standard Big Bang phase of the Universe. At high 
frequencies, however, the photon's dispersion relation re- 
flects the structure and properties of the crystal. The 
crystal breaks Lorentz invariance, and like in the case 
of phonons, the photon's speed becomes frequency de- 
pendent. It is in this phase that the horizon and flatness 
problems are solved. The fact that the Universe is pushed 
towards a colder phase, for which non-commutativity and 
violations of Lorentz invariance become negligible, can be 
seen as an explanation for the emergence of the contin- 
uum and symmetries we perceive today. 

In a future publication [23] we shall go further and 
prove that thermal fluctuations in deformed radiation 
may translate into a Harrison- Zeldovich spectrum of ini- 
tial conditions for structure formation. Crucially this 
rules out all but a well-dcflned class of quantum defor- 
mations. This important result provides a thermal al- 
ternative to the usual quantum origin of the structure 
of the Universe (which may or may not have its prob- 
lems [30,31]). It also, for the first time, allows VSL to 
become a proper model of structure formation. Another 
avenue to be studied further concerns the cosmological 
constant problem. Moffat has recently [25] shown how 
a resolution may emerge from non-local field theories. 
Non-commutativity imposes just such a non-locality. 

In a somewhat different direction in [32] we describe 
our findings for model 3, with 7 > 1, a range left un- 
touched in this paper. We find that ordinary thermal 
radiation subject to this type of non-commutativity may 
drive inflation. We thus realize the inflationary scenario 
without the aid of an inflaton field. As the radiation 
cools down below Planck's temperature, infiation grace- 
fully exits into a standard Big Bang universe, dispensing 
with a period of reheating. Curiously in this regime there 
is no VSL, as the color temperature is found to saturate. 

Recently, progress in String/M theory and Spin Net- 
works have realized a non-commutative phase. We 
have conjectiircd a non-commutative version of the 
Friedmann- Walker space-time in our realization of the 
equivalence of VSL and k deformed non-commutativity 
in order study the dynamics the early universe. We ex- 
pect that this type of non-commutativity to come from 
a non-constant B field in a curved background in String- 
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theory. Wc appreciate that this is beyond the state of 
the art in string-theory. However, it would be impor- 
tant to find a more concrete realization of NCFRW from 
string/M theory. 
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